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The charge dynamics of the copper oxide materials in the underdoped and optimal 
doped regimes is studied within the framework of the fermion-spin theory. The 
conductivity spectrum shows the non-Drude behavior at low energies and unusual 
midinfrared peak, and the resistivity exhibits a linear behavior in the temperature, 
which are consistent with experiments and numerical simulations. 
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Since the discovery of the copper oxide superconductors, a significant body of reliable 
and reproducible data has been accumulated by using many probes [|1] , which show that the 
most remarkable expression of nonconventional physics of copper oxide materials is found 
in the normal-state. The normal-state properties exhibit a number of anomalous properties 
in sense that they do not fit in the conventional Fermi-liquid theory, and some proper- 
ties mainly depend on the extend of dopings 0. The undoped copper oxide materials are 
antiferromagnetic Mott insulators, and the antiferromagnetic long-range-order (AFLRO) 
disappears when the hole doping concentration exceeds some critical value (about 5%) [||. 
The charge response in the underdoped and optimal doped regimes, as manifested by the op- 
tical conductivity and resistivity, has been extensively studied experimentally in the copper 
oxide materials, as well as theoretically within some strongly correlated models • From 
experiments testing normal-state properties, the optical conductivity shows the non-Drude 
behavior at low energies and anomalous midinfrared band in the charge-transfer gap, the 
resistivity exhibits a linear behavior in the temperature in the optimal doped regime and 
a nearly temperature linear dependence with deviations at low temperatures in the under- 
doped regime 000] • These optical properties indicate that the strong correlations are very 
important to the electronic structure. 

Among the microscopic models the most helpful for the discussion of the normal-state 
properties in the copper oxide materials is the t-J model acting on the space with no doubly 
occupied sites ||. In order to account for real experiments under the t-J model, the crucial 
requirement is to impose the electron local constraint ||. Recently a fermion-spin theory 
based on the charge-spin separation is proposed |L0| to incorporated this constraint. In this 



approach, the electron on-site local constraint for single occupancy is satisfied even in the 
mean-field approximation (MFA). Within the fermion-spin theory, it has been shown [DJ 
that AFLRO vanishes around doping 5 = 5% for the reasonable value of the parameter t/J = 
5, and the mean- field theory |12j in the underdoped and optimal doped regimes without 



AFLRO has been also developed, which has been applied to study the photoemission and 
electron dispersion. In this paper, we try to study the charge dynamics of the copper oxide 



materials within this fermion-spin theory. We follow the fermion-spin formulism HlO| , |T2| and 
decompose the electron operators as = h\S^ and = h\Sf, where the spinless fermion 
operator hi keeps track of the charge (holon), while the pseudospin operator Si keeps track 
of the spin (spinon). The mean-field order parameters are defined |L^] as \ — (SfST+n) = 

(&i Si+7)}> Xz = (SiSi+ v )i C = (V^ 2 ) ^rj,v' (&T+ rj^i+n' ) > @z = (V^ 2 ) ^T),if (^i+n^i+rf) ' anC ^ 

= (h\hi +r) ) with f) = ±x, ±y, and Z is the number of nearest neighbor sites. In this case, 
the low-energy behavior can be described [E| by the effective Hamiltonian H = H t + Hj 



with 

H t = -t y: hrt +v (s+sr +v + srs+ +v ) + h. c . + E h\h h (i a ) 

Hj = Jeff Y^2^^ +r > + + ( lb ) 

where J e ff = J[(l — <5) 2 — 2 ], and /x is the chemical potential which enforce (h\hi) = 5. 

The spinon and holon may be separated at mean-field level, but they are strongly 
coupled beyond MFA due to fluctuations |T3]| . For discussing the optical conductivity 
and resistivity, we now need to consider the dynamical effects beyond MFA. Since an 
electron is represented by the product of a holon and a spinon in the fermion-spin the- 
ory [ID,r^|, the external field can only be coupled to one of them. According to Ioffe- 



Larkin combination rule fl4[], we can separately calculate the contributions to the con- 
ductivity from holons (Th{u) and spinons <7 s (u), and the total conductivity is given by 
(j _1 (cj) = o"^ 1 (cj) + cr" 1 ^). In the present theoretical framework, the current densities of 
spinons and holons are defined as j s = te<p ^2 ir) f)(Sf S^ v +S~ S^_ v ) and jh = 2te% vh\ +r] hi, 
respectively. A formal calculation for the spinon part shows that there is no the direct 
contribution to the current- current correlation from spinons, but the strongly correlation 
between holons and spinons is considered through the spinon's order parameters x j Xz 
, C, and C z entering in the holon part of the contribution to the current- current corre- 
lation, which means that the holon moves in the background of spinons, and the cloud 
of distorted spinon background is to follow holons, therefore the dressing of the holon by 
spinon excitations is the key ingredient in the explanation of the transport property. In 



this case, the conductivity is expressed as cr(u) = — lmllh(uj)/uj with the holon-holon cor- 
relation function Ii h {iu n ) = -(2texZ) 2 (l/N) £ fc lsk( l /P) £iu; m , 9(k, iu m > + iu n )g(k,iu m t), 
where jsk = (l/2)(sin/c x + smk y ) and the holon Green's function is defined as gij(r — r') = 
— (T T hi(r)hj(r')) . In this paper, we limit the spinon part to the first-order (mean-field level) 
since some physical properties can be well described at this level []12||, and there is no di- 



rect contribution to the charge dynamics from spinons as mentioned above. However, the 
second-order correction for the holon is necessary for a proper description of holon motion 
due to antiferromagnetic fluctuations. The mean-field spinon Green's functions D(°\k,uj) 
and D^(k,uj) and mean-field holon Green's function g(°\k,uj) have been given in Ref. ||12|| . 
The second-order holon self-energy diagram from the spinon pair bubble is shown in Fig. 1. 
Since the spinon operators obey the Pauli algebra, we map the spinon operator into the CP 1 
fermion representation or spinless-fermion representation in terms of the 2D Jordan- Wigner 



transformation [15| for the formal many particle perturbation expansion. After then the 
spinon Green's function in the holon self-energy diagram shown in Fig. 1 is replaced by the 
mean- field spin Green's function D^(k, uj), and therefore the second-order holon self-energy 
is evaluated as 

1 



(k,uj) = {Ztf— J2(%'-k + lp'+ P +k) 2 B pl B p+pl x 



, "f(^)[»bK') - n B {up + p')] - n B {u p+p i)n B (-u p ,) 
uj + uj p + p i — ujp> — £ p+ fc + i0 + 
n F (£,p+k)[n B {uj p+p >) - n B (-uJp')] + n B {uj p <)n B {uj p+p > 



uj + L0 pl + u p+p/ - £ p+fe + i0 + 
n F (^ p+k )[n B (uj p+p ') -n B (-u)p>)} - n B (-uj p >)n B (-uj p + p/ ) \ 
uj — uj p+p i — ujpi — £ p+ fc + i0 + J 

where 7fc = (1/Z) £„e^, e = 1 + 2t<j>/J eff , B k = ZJ eff [{2e Xz + x)lk ~ (e X + 2 Xz )]/uj k , 
n F{£,k) and n B (uJk) are the Fermi and Bose distribution functions, respectively, the mean- 
field holon excitation spectrum = 2Z\t r )k + and the mean-field spinon excitation 
spectrum u k is given in Ref. ||12| . Then the full holon Green's function is obtained as 
g^ 1 {k,uo) = g^^ 1 {k,uo) — £^ (k,u). We emphasize that the local constraint of the t-J 
model has been treated exactly in the mean-field theory [O, and it is still satisfied in the 



above perturbation expansion based on this mean-field theory. 

With the help of the full holon Green's function, the optical conductivity can be obtained 

as 

/ \ r?\2 1 2 f°° du)' . . \ a ft n n F (uo' + uo) - n F (uo') 

with the holon spectral function Ah(k,u) = —2lmg(k,uo). We have performed a numerical 
calculation for this optical conductivity (3), and the results at the doping 5 = 0.06 (solid 
line) and 5 = 0.12 (dashed line) for the parameter t/J — 2.5 with the temperature T = 
are shown in Fig. 2, where the charge e has been set as the unit. For comparison, the 
result from numerical simulation based on the Lanczos diagonalization || of a 4 x 4 cluster 
with a single doped hole (5 ~ 0.06) for t/J — 2.0 is also shown in Fig .2 (dot-dashed line). 
Although the infrared properties of the copper oxide materials are very complicated, some 
qualitative features such as away from half-filling, weight appears inside the charge-transfer 
gap of the undoped materials, defining the midinfrared band and the conductivity decays 
as — > 1/uj at low energies, seem common to all copper oxide materials [0,0,0. The present 
theoretical result shows that there are a low-energy peak at u < 0.5t separated by a gap or 
pseudogap ~ 0.5t from the broad absorption band (midinfrared band) in the conductivity 
spectrum. The analysis of the result indicates that the conductivity decays is — > 1/u (non- 
Drude falloff) at low energies (u < 0.5t), the midinfrared peak is doping dependent and the 
peak is shifted to lower energy with increased doping, which is in qualitative agreement with 
the experiment 0] and numerical simulation |5]||. 

It has been shown that the Drude weight can be used as an order parameter for metal- 



insulator transitions [16], which converges exponentially to zero for the insulator, and 



nonzero constant for the metal. Following the standard discussions |p.7|| , the Drude weight 
D is obtained in the present case as 

r-, ,. /r , , n2 1 • 2, f°° duo duo' n F {u') - n F (uo) . A . 

D n = ^—i- + ,hm (2 X t ^E sin ^ / ^= 1 r-A h (k,uo)A h (k,uo') \ ' \ J , (4) 

4 tuj n ^0 iV £ J-oc 2,71 2n iuo n + u — uo 

where D n = D/27re 2 , and the total kinetic energy per site (T) = 8tx4>- The numerical 
results of eq. (4) for t/J — 2.5 (solid line) is plotted in Fig .3. The result shows that the 



Drude weight grows rapidly with increasing doping 5, and the doping dependence is nearly 
linear in the underdoped and optimal doped regimes, which is qualitative consistent with 
the exact diagonalization result |18[ . 

The resistivity is expressed as p = 1/crdc, where the dc conductivity o<ic is obtained from 
Eq. (3) as o^c — h m t^o (?(uS). This resistivity has been evaluated numerically and the results 
in the doping 5 = 0.06 and 5 = 0.10 for t/J— 2.5 are plotted in Fig. 4(a) and Fig. 4(b), 
respectively, which show that in the underdoping, the resistivity indeed exhibits a nearly 
temperature linear dependence with deviations at low temperature, and a very good linear 
behavior in the optimal doping, which also is in agreement with the experiments fl]|||4]]. 

In the present theory, the basic low-energy excitations are holons and spinons. However, 
our results show that the unusual non-Drude type optical behavior at low energies and 
midinfrared band, and linear temperature dependence of the resistivity are mainly caused 
by the holons in the copper oxide sheets, which are strongly renormalized because of the 
interactions with fluctuations of the surrounding spinon excitations. The 1/uj decay of the 
conductivity at low energies is closely related with the linear temperature resistivity, this 
reflects an anomalous frequency dependent scattering rate proportional to u instead u 2 as 
would be expected for the conventional Fermi-liquid theory, which is consistent with the 
Luttinger-liquid theory and phenomenological marginal Fermi-liquid theory |20| Our 
study seems to confirm that the strongly correlated t-J model is possibly sufficient for the 
understanding of the normal-state properties of the copper oxide materials. 

In summary, we have studied the charge dynamics of the copper oxide materials in the 
underdoped and optimal doped regimes within the t-J model. It is shown that the strongly 
correlated renormalization effects of the holon motion due to spinon fluctuations are very 
important for the charge dynamics. The optical conductivity, Drude weight and resistivity 
are discussed, and the results are qualitative consistent with the experiments and numerical 
simulations. 
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FIGURES 

FIG. 1. The holon's second-order self-energy diagram. The solid and dashed lines correspond 
to the holon and spinon propagators, respectively. 

FIG. 2. The optical conductivity at the doping 5 = 0.06 (solid line) and 5 = 0.12 (dashed line) 
for t/J = 2.5 with the temperature T = 0. The dot-dashed line is the result from the numerical 
simulation on a 4 x 4 cluster with a single doped hole for t/J = 2.0 [||. 

FIG. 3. Drude weight, in the unit of t, as a function of doping 6 for t/J = 2.5. 

FIG. 4. The electron resistivity at the parameter t/J = 2.5 for (a) the doping 5 = 0.06 and (b) 
5 = 0.10. The small wiggles are the finite-size effect. 
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